Ab-initio study on the possible doping strategies for M0S2 monolayers 
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Density functional theory is used to systematically study the electronic and magnetic properties 
of doped M0S2 monolayers, where the dopants are incorporated both via S/Mo substitution or as 
adsorbates. Among the possible substitutional dopants at the Mo site, Nb is identified as suitable 
p-type dopant, while Re is the donor with the lowest activation energy. When dopants are simply 
adsorbed on a monolayer we find that alkali metals shift the Fermi energy into the M0S2 conduc- 
tion band, making the system n-type. Finally, the adsorption of charged molecules is considered, 
mimicking an ionic liquid environment. We find that molecules adsorption can lead to both n- and 
p-type conductivity, depending on the charge polarity of the adsorbed species. 

PACS numbers: 



I. INTRODUCTION 

In recent years two dimensional (2D) materials have 
attracted a growing interest due to their potential for 
future nano-electronics applications, owing to their un- 
usual physical, optical and electrical properties arising 
from the quantum confinement associated to their ul- 
tra thin structure^. These considerations apply partic- 
ularly well to layered transition metal di-chalcogenides 
(LTMDs) in which a multitude of electronic states have 
been observed, namely metallicity, semi-conductivity and 
charge density waves. Given the richness in the electronic 
properties it is not a surprise that this materials class has 
been called out as the ideal platform for a multitude of 
application^213] Pqj. example, the prototypical LTMD, 
molybdenum di-sulfide M0S2, has been widely explored 
as lubricantpl, catalysiP and as lithium ion battery an- 
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Structurally M0S2 consists of covalently bonded S-Mo- 
S 2D hexagonal planes (monolayers), which in the bulk 
are bound together in a layered structure by weak van 
der Waals forces . Interestingly, the electronic proper- 
ties of bulk M0S2 show a strong dependence on the layer 
thickness^. Single M0S2 monolayers are particularly in- 
triguing. These display a direct band-gap of 1.9 eV, 
which makes thenr a senriconducting alternative to gap- 
less graphene. Indeed the band-gap of graphene can be 
opened by fabricating nanoribbons- or by depositing it 
over a suitable substrate^'^. However this comes at the 
cost of deteriorating the carrier mobility due to edge and 
impurity scattering^. In contrast, the absence of dan- 
gling bonds, the high crystallinity, and the low dimen- 
sionality, make the performance of LTMDs comparable 
to those of currently existing Si transistors at the scaling 
limilJiSHlll. M0S2 monolayer based transistors have been 
recently demonstrated to operate at room temperature, 
with a mobility of at least 200 cm^/Vs and on/off cur- 
rent ratios of 10* with low standby power dissip ation^. 
Interestingly, both n-type^^HiSl and p-typ^iSMl conduc- 
tivities have been reported for ultra-thin M0S2 layers, 
depending on the experimental details. The origin of 
such diverse conducting properties remains to date far 



from being clear. 

The possible creation of Mo and/or S vacancies in 
M0S2 monolayers during the growth cannot be used as 
a doping strategy, since those vacancies always produce 
deep trap states in the M0S2 monolayer band-gap22l. In a 
recent study we have investigated how unintentional dop- 
ing of the substrate holding a M0S2 monolayer may deter- 
mine its conductivity2lI. Here we explore possible strate- 
gies for doping directly into the M0S2 monolayer. Al- 
though there are studies available on specific dopants for 
single layer MoSj^^HSU^ ^q q^j. knowledge there is not yet a 
systematic investigation comparing the doping properties 
of the various defects, when they are incorporated in ei- 
ther a substitutional position or they are adsorbed on the 
surface. Moreover, it has been reported that in a liquid- 
gated electric double layer transistor, the ions can modu- 
late the electronic properties from insulating to metallic 
through electrostatic induced carrier doping in the trans- 
port channel at a finite gate bias^H. Compared with other 
dopants, such as metals and light atoms, ionic liquids 
have several advantages as surface dopants: a large vari- 
ety of such ionic complexes are available, they can accu- 
mulate more carriers than conventional solid state gated 
transistors, and they do not induce structural disorder 
like substitutional dopants. In fact, ambipolar transistor 
operation has been achieved using thin flakes of M0S2 in 
an ionic liquid gated environment^^''. The present work 
explores also this possibility. 

Here we systematically study the effects of dopants on 
the electronic structure of a M0S2 monolayer, by calcu- 
lating the formation energies and the electronic proper- 
ties of halogens, non-metals, transition metals and alkali 
metals added in various geometrical conflgurations. Our 
aim is that of identifying potential candidates for mak- 
ing M0S2 monolayers either n- or p-type in a controlled 
way. Then we present an investigation of the electronic 
structure of M0S2 monolayers when two molecular ions, 
NH4" and BF4 , are adsorbed at the surface. The paper is 
orgarrized as follows. After a brief discussion of the com- 
putational methods used we turn our attention to the 
case of substitutional doping by looking at both at the 
S and the Mo site. Next we consider the case of doping 



by adsorption of alkali and molecular ions. At the end 
of the results section we briefly discuss the robustness of 
our calculations against the particular choice of exchange 
and correlation functional and finally we conclude. 



II. METHODOLOGY 

In order to investigate the electronic properties of a 
M0S2 monolayer doped with impurities, ab-initio cal- 
culations are performed using density functional the- 
ory2ZB3 within the local spin density approximation 
(DFT-LSDA) for the exchange and correlation poten- 
tial. In particular we consider the Ceperly- Alder LSDA 
parametrizatiorPSlas implemented in the SIESTA code^Sl. 
In our calculations double-^ polarizecpi' numerical atomic 
orbitals basis sets are used for all atoms, and the 
TrouUier-Martins scheme is used for constructing norm- 
conserving pseudopotentials^. A 5x5 hexagonal super- 
cell [see Fig. [Ha)] with (15.66x15.66) A^ lateral dimen- 
sions is constructed and doping is introduced by replac- 
ing/adding a single atom in the supercell. This corre- 
sponds to simulating a Moi_a;Aa;S2-j;Bj, periodic system 
with X — 4% {y — 2%) when the doping is at the Mo 
(S) site. An equivalent plane wave cutoff of 250 Ry is 
used for the real space mesh and the Brillouin zone is 
sampled over a 5x5x1 Monkhorst-Pack fc-grid. Periodic 
boundary conditions are applied and a vacuum layer of at 
least 15 A is placed above the monolayer to minimize the 
interaction between the adjacent periodic images. A tem- 
perature of 300 K is used when populating the electronic 
states with a Fermi distribution. The relaxed geometries 
are obtained by conjugate gradient, where all the atoms 
in the supercell are allowed to relax until any force is 
smaller than 0.02 eV/A. 

In order to verify that the calculated impurity level 
alignments are robust against the choice of exchange cor- 
relation functional, we repeat the calculations for the 
main results by using the screened hybrid functional of 
Heyd-Scuseria-Ernzerhof (HSEOe^P. All DFT calcula- 
tions based on the IISE06 functional are carried out 
with the projector augmented wave (PAW) pseudo poten- 
tial plane-wave method^ as implemented in the VASP 
code^^. A 3x3x1 Monkhorst-PaclP^ fc-point grid is em- 
ployed and the plane wave energy cutoff is 500 eV. 

The formation energy of a particular substitutional 
dopant, inform, is defined as 



£^form = £;tot(MoS; + X) + £;bulk(h0St) 
~Stot(M0S2)-^bulk(X), 



(1) 



where £'tot(MoS2-l-X) is the total energy of the system 
including the substitutional atom X, £'tot(MoS2) is the 
total energy of the corresponding pristine M0S2 mono- 
layer, while £^buik(X) and ii^buik(host) are respectively the 
energies of the substitutional atom X and of the substi- 
tuted Mo (S) host atom in their bulk (diatomic molecule) 
forms. In contrast the formation energy for an adsorbate. 




FIG. 1: (Color online) Structure of the supercell used for the 
calculations. In (a) we show a 5x5 M0S2 monolayer supercell 
with the possible substitutional dopants. Mo is substituted by 
a transition metal (green dot): Y, Zr, Nb, Re, Rh, Ru, Pd, Ag, 
Cd; S is substituted by nonmetals and halogens (red dot): P, 
N, As; F, CI, Br, I. In (b) we display the same supercell where 
we indicate the possible adsorption sites defined in the text. 
Light gray spheres indicate Mo atoms, light yellow spheres 
indicate S. 



-Eadsj can be written as 

£;ads = ^tot(MoS2 + Y) - £;tot(M0S2) - Sbulk(Y) , (2) 

where i?tot(MoS2-l-Y) refers to the total energy when the 
adsorbate Y is attached to M0S2, and £'buik(Y) is the 
energy of an adsorbate Y in its bulk form. In order to find 
the most stable configuration for adsorption, we consider 
four possible positions labeled as follows: Tg (adsorbate 
on top of S), Tmo (adsorbate on top of Mo), A (adsorbate 
above the center of the hexagonal ring of M0S2), and 
B (adsorbate above the middle of the Mo-S bond) [see 
Fig.Hb)]. 



III. RESULTS AND DISCUSSION 

Before discussing the formation energy and the elec- 
tronic structure of the possible dopants let us here briefly 
review the electronic properties of a single M0S2 mono- 
layer. As the M0S2 thickness is decreased from bulk 
to a few layers, the valence band minimum shifts from 
half way along F-K line towards K ^. For a mono- 
layer the band-gap becomes direct at K (see Fig. [2|, a 
transition which has been recently observed experimen- 
talljSI. The computed LSDA band-gap of 1.86 eV is in 
good agreement with the experimental optical band-gap 
of 1.90 e\^, although such an agreement has to be con- 
sidered fortuitous. In fact the absorption edge of an op- 
tical excitation measures the energy difference between 
the quasi-particle band-gap and the exciton binding en- 
ergy. In M0S2 monolayers this latter is of the order of 
1 eV, as confirmed recently by many-body calculations^Zl. 
Thus, one expects that the true quasi-particle spectrum 
has a band-gap of approximately 2.9 eV, in good agree- 
ment with that computed with the GW scheme, either 
at the first order leveP^ (2.82 eV) or self-consistentljEH 
(2.76 eV). At the end of this section we will discuss how 
band-gap corrections, via the HSE06 functional, affect 
our main results. In any case, for both bulk and mono- 
layer the band-structure around the band-gap is derived 
mainly from Mo-Ad orbitals (see the orbital projected 



> ' 



PDOS (states/eV/spin) 






FIG. 2: (Color online) Band-structure (left) and DOS (right) 
of a pristine (undoped) M0S2 monolayer. The red curve indi- 
cates the DOS projected on the Mo atoms. 



Impurity 


m(^s /impurity) 






Eform(eV) 


F, CI, Br 




1.00 




0.60 


(F), 


1.93 (CI), 2.16 (Br) 


I 




0.00 








2.53 


N, P 




1.00 






2.90 


(N), 1.89 (P) 


As 




0.00 








1.91 


Re 




1.00 








2.05 


Ru 




2.00 








3.05 


Rh 




3.00 








4.15 


Pd 




4.00 








5.63 


Ag- 




2.57 








7.28 


ed 




1.58 








6.11 


Nb, Zr, Y 




0.00 




-0.19 


(Nb) 


-0.48 (Zr), 0.52 (Y) 



TABLE I: The theoretically calculated magnetic moment, m, 
and the formation energy, -Eform of different substitutional 
dopants in the M0S2 monolayer. 



density of state, PDOS, in Fig. [2]), although there are 
also small contributions in the valence band from the S- 
3p ones. 



A. Substitutional doping 

Substitution at the S site 

We begin our analysis by substituting a surface S atom 
with elements taken from the halogen family, namely F, 
CI, Br and I. These are expected to act as a source of 
n-type doping for M0S2, since they have one additional 
p electron with respect to S. The DOSs for the supercell 
including one atom of the halogen family are presented in 
Fig. [3j As a representative system, we discuss in details 
results for Cl-doping, since the other halogens present 
a similar electronic structure. A Cl-doped M0S2 mono- 
layer has a magnetic ground state, where an occupied 
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FIG. 3: (Color online) DOSs for a 5x5 M0S2 supercell in 
which one S atom is replaced by (a) F, (b) CI, (c), Br and (d) 
I. Negative DOS values refer to minority spins, while positive 
are for the majority. The blue dashed line indicates the Fermi 
energy, while the colored shaded areas indicate the DOS pro- 
jected over the dopants. Note that all the DOSs are aligned 
to have a common Fermi level, Ef = 0. 



defect level is formed at about 0.4 eV below the conduc- 
tion band minimum (CBM). The corresponding minority 
state is located above the CBM and is empty. These de- 
fect levels originate from the hybridization between the 
Cl-3p and the Mo-Ad states. Similarly to the CI substi- 
tution, for both F and Br dopants (isoelectronic to CI) 
the system is para-magnetic, having a magnetic moment 
of approximately 1 fj-B- In the case of F substituting 
S, the spin-splitting is larger than that of CI, whereas 
it is smaller for the Br substitutional case. Eventually, 
when S is replaced by iodine, the impurity becomes non- 
magnetic, although the impurity level is still located well 
below (0.3 eV) the CBM [see Fig.jSJJd)]. Note that hav- 
ing an impurity presenting a ground state with finite 
magnetic moment has nothing to do with diluted ferro- 
magnetism in p-type systems, as sometimes erroneously 
claimed in literature. This simply indicates that the ad- 
ditional electron remains unpaired and localized around 
the impurity'^^. The dependence of the spin-splitting of 
the impurity level on the atomic number simply follows 
the well-known trend for the exchange-correlation inte- 
gral"*". 

We then investigate the possibility of p-type doping 
obtained by replacing a S atom with a group V element 
of the periodic table, namely N, P and As (see Fig. 4]). In 
the case of N the supercell has a magnetic ground state 
[see Fig. Ka)], with a magnetic moment of IfiB and a 
spin-splitting of about 0.20 eV. In both the spin chan- 
nels gap states are introduced above the valence band 
maximum (VBM) , and these are formed by hybridization 
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FIG. 4: (Color online) DOSs for a 5x5 M0S2 supercell in 
which one S atom is replaced by (a) N, (b), P and (c) As. 
Negative values refer to minority spins, positive values to ma- 
jority spins. The blue dashed line indicates the Fermi energy, 
the colored shaded areas indicate the DOS projected on the 
dopants. Note that all the DOSs are aligned to have a com- 
mon Fermi level, Ey = 0. 



between the N-2]3 and the Mo-4d valence orbitals. Sim- 
ilarly, when substituting S with P, the system is again 
magnetic, but the defect states are rigidly shifted closer 
to the VBM of the pristine M0S2 monolayer in both the 
spin channels. The spin-spUtting remains similar to that 
of N. Finally for the As substitutional case, the defect 
states shift even more towards the VBM. Additionally, 
the spin-splitting vanishes and the system becomes a p- 
type semiconductor. The partially unoccupied impurity 
band however shows very little dispersion and it is sep- 
arated from the valence band by about 0.08 eV. Such 
an impurity band is formed mainly by As-4p and Mo-4d 
orbitals. 

To summarize our results for S substituting impuri- 
ties, we find that most dopants create localized, spin- 
polarized, gap states. When increasing the dopant 
atomic number the spin-splitting reduces and the states 
move towards the VBM/CBM. Among the many dopants 
investigated As appears to give the most promising DOS 
for p-type doping, since the associated diamagnetic im- 
purity band is very close to the VBM. In the case of 
n-type doping, all substituent dopants have donor states 
located rather far from the M0S2 CBM and therefore 
they are rather localized. In addition to these consid- 
erations a look at Table |l) where we present the forma- 
tion energy of the various dopants, reveals that all the S 
substituents possess a rather large formation energy, i.e. 
they are unlikely to form under equilibrium thermody- 
namical conditions. The only viable formation channel 
is offered by filling S vacancies, which have been recently 
demonstrated to form with relatively ease^. Therefore, 
it appears that a strategy for doping at the S site may 
be that of growing S poor samples and then of filling 
the vacancies with an appropriate donor/acceptor sulfur 
replacement. 



Substitution at the Mo site 

Next we consider substitutional doping at the Mo site 
with different transition metal atoms having a growing 
number of electrons in the d shell. Specifically we con- 
sider all the Ad elements in the periodic table ranging 
from Y to Cd with the only exception of radioactive Tc, 
which is replaced by Re (5d). In general we find that for 
all the substitutional atoms having a d occupancy larger 
than Mo, the ground state is spin-polarized, while Y, Zr 
and Nb have a diamagnetic ground state (see Table IT]). 

Let us again discuss the case of n-type doping first. 
When Re replaces Mo, one extra electron is added to the 
supercell. This remains unpaired and occupies a majority 
spins gap state [Fig.lsla)], located 0.3 eV below the CBM, 
with the Fermi energy, E-p located between such state and 
the CBM. As a consequence the supercell has a magnetic 
moment of 1 /is and the DOS for this situation resem- 
bles that of CI substituting S [Fig. Islb)]. If one forces 
a non-spin-polarized solution for Re-doped M0S2 mono- 
layer, the donor level is created at about 0.2 eV below 
the CBM. However, the magnetic ground state is lower in 
energy than the non-magnetic one by 82 meV/supercell, 
i.e. a diamagnetic solution is not stable. Note that our 
non-spin-polarized result is consistent to that reported 
for similar a calculatiorP^l, where a Re donor level at 
0.19 eV below the CBM suggested that Re could be used 
as n-type dopant in M0S2 nanotubes. Such result thus 
seems to be robust against spin polarization (the spin po- 
larized impurity state is at 0.3 eV below the VBM). Note 
that single crystals Re-doped M0S2 have been grown in 
the past by chemical vapor deposition -. 

When one then looks at other transition metal dopants 
having a d-orbitals occupancy larger than Re, an increas- 
ing number of gap states is formed. These are progres- 
sively occupied so that the excess electrons do not spill 
into the M0S2 conduction band (see Fig. [5]). As a conse- 
quence the Fermi energy always lies in the M0S2 band- 
gap, and in fact it moves deeper into the band-gap as 
the d shell filling of the dopant increases. The excess 
electrons first fill up the majority spin states of the host 
d-orbitals, until the magnetic moment of the impurity 
reaches the largest value of 4 fiB for Pd. Then the mag- 
netic moment decreases for Ag and Cd as the additional 
electrons start to populate the minority spins states (see 
Table [!]). 

Next we move to study the possibility of obtaining p- 
type doping by replacing Mo with Nb, Zr and Y, which 
have respectively one, two and three electrons less than 
Mo. The case of Nb seems to be particularly attrac- 
tive. The inclusion of Nb into M0S2 changes little the 
bond lengths and also the main DOS. The only notable 
effect is the shift of the Fermi energy below the VBM 
due to the one electron removal [see Fig. |6][a)]. The 
newly created defect states are rather delocalized and 
the charge excess spreads out up to the third nearest 
neighbor Mo atoms. The states around the VBM mainly 
originate from hybridized d orbitals of Nb and Mo. The 
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FIG. 5: (Color online) DOSs for a 5x5 M0S2 supercell in 
which one Mo atom is replaced by (a) Re, (b) Ru, (c) Rh, (d) 
Pd, (e) Ag and (f) Cd. The blue dashed line marks the Fermi 
energy, the colored shaded areas indicate the DOS projected 
onto the dopants. Note that all the DOSs are aligned to have 
a common Fermi level, _Bf — 0. 




FIG. 6: (Color online) DOSs for a 5x5 M0S2 supercell in 
which one Mo atom is replaced by (a) Nb, (b) Zr and (c) 
Y. The blue dashed line indicates the Fermi energy, the col- 
ored shaded areas indicate the DOS projected on the dopants. 
Note that all the DOSs are aligned to have a common Fermi 
level, Ef = 0. 



Atom-Atom 5E (moV) 



Nb-Nb 
Cs-Cs 
Nb-Cs 



-157 
166 
-218 



TABLE II: Pairing energy for two impurities doping a M0S2 
monolayer. Calculations are based on total energy difference 
between a nearest neighbor and a separated geometry in a 
6x6 supercell. Negative 5E indicate a tendency to clustering. 



valence band now looks sufficiently dispersive in our DFT 
band-structure, and therefore the mobility is expected to 
be rather large. Importantly the formation energy for 
Nb-doping is -0.19 eV, in good agreement with the for- 
mation energy of -0.21 eV^ll, obtained theoretically for 
Nb in bulk M0S2. Overall our results suggest that Nb 
may be a promising candidate as p-type dopant in M0S2 
monolayers. Note that an experimental study^^ shows 
that Nb-substituted (concentration range 15-25 %) M0S2 
nanoparticles can be synthesized, and that they also ex- 
hibit p-type character. Finally, for all the other two p- 
type dopants considered (Y and Zr) the ground state is 
also non- magnetic (see Table ll|. However, in contrast to 
Nb doping, the defect state becomes less hybridized with 
the VBM and produces split off acceptor-levels above the 
VBM for both Zr and Y [see Fig.jejb) and (c)]. 

We conclude this section by evaluating the pairing en- 
ergy between two Nb atoms placed in a 6x6 supercell. In 
Table [TT| we list the energy difference SE between the con- 
figuration where the two impurities are placed at nearest 
neighbor positions and when they are placed as far as 
possible in the supercell. The result suggests that in case 



of Nb pairs the dopants tend to stay close to each other 
in the M0S2 monolayer {SE < 0). Thus, based solely 
on the pairing energy, Nb dopants would form clusters. 
Clustering however is likely to be inhibited in by the large 
energy barrier for Nb diffusion in the M0S2 plane. 

Summarizing the situation for transition metal doping 
at the Mo site we find that, when the dopant has more d 
electrons than Mo, donor states are created deep inside 
the M0S2 band-gap (at least 0.2 eV below the CBM), 
with Re being the dopant with the smallest activation 
energy. However, in all cases the formation energies are 
large (Table |l| . On the contrary, p-type doping obtained 
by substituting Mo with transition metals such as Nb 
and Zr creates acceptor states just at the VBM, and the 
formation energies are also small (Table IT]). In fact, for 
both Nb and Zr, inform is negative, indicating that substi- 
tutional doping will form spontaneously in M0S2 mono- 
layers (note that this is relative to the bulk reference for 
the dopant). Intriguingiy, also for the case of transition 
metal doping at the Mo site the possibility of filling va- 
cancies remains open. In fact a recent experimenipSI has 
provided evidence for the formation of Mo vacancies in 
bulk M0S2 via proton irradiation. 




FIG. 7: (Color online) DOSs for a 5x5 M0S2 supercell in 
which an atom chosen between (a) Cs, (b) K, (c) Li and (d) 
H is adsorbed on the M0S2 surface. The blue dashed line 
indicates the Fermi energy, the colored shaded areas indicate 
the DOS projected on the adsorbates. Note that all the DOSs 
are aligned to have a common Fermi level, Ep — 0. 



Adatom E 


.ds (eV) 


Cs 


-0.79 


K 


-0.82 


Li 


-0.98 



TABLE IIL Theoretically calculated adsorption energy for 
different alkali metal adsorbed on the M0S2 monolayer. 



B. Doping by adsorption 

The results of the previous section show that neither 
Mo substitution with transition metals nor S substitu- 
tion with non-metal elements appear as promising strate- 
gies for obtaining shallow donor states (p-doping is much 
more promising with Nb). Another possible route for 
obtaining n-type M0S2 monolayers is by adsorbing H or 
alkali metals such as Li, K and Cs. This possibility is 
explored here. Then, in the second part of this section 
we consider adsorbed molecular ions as potential dopants 
with both n- and p-type character. 



Alkali atom adsorption 

Among all the possible adsorption sites for alkali met- 
als on the M0S2 surface, we find that the Tmo one [see 
Fig. WCb)] is the most energetically favorable. This is 
also suggested experimentalljSZl as well as predicted by a 
previous theoretical calculationa^Sl. Most importantly all 
the adsorption energies are found large and negative (see 
Table III I indicating thermodynamical stability. 

We start by considering Cs adsorption. When Cs is 
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FIG. 8: (Color online) DOS for a M0S2 supercell, where a 
single Mo is replaced by Nb, and where a Cs atom is adsorbed. 
The blue dashed line indicates the Fermi energy, Ep — 0, the 
colored curves indicate the DOS projected on the dopants. 



adsorbed on a M0S2 monolayer there is no significant 
change in the geometry, i.e. the ion binds without dis- 
torting the lattice of the host. The DOS presented in 
Fig.lWa) shows that the available Cs 6s electron is trans- 
ferred to the CBM of the M0S2 monolayer, occupying the 
Mo-Ad orbitals. In this case there is no split off at the bot- 
tom of the conduction band and the DOS projected over 
the Cs 6s orbitals appears uniformly spread over a 1.5 eV 
energy window. This essentially means that Cs adsorbed 
over a M0S2 monolayer acts as perfect donor. A very 
similar situation is found for both K and Li, with only 
minor quantitative differences in the spread of the DOS 
associated with the s shell of the dopant [see FigJTJ^b) 
and FigJTFc)]. In contrast, the adsorption of H produces 
a spin-split state 1 eV below the CBM, i.e. roughly at 
midgap of the LSDA band-gap [see FigJTl. As such we 
conclude that H is not a shallow donor for M0S2, in good 
agreement with the recent theoretical results^. 

The pairing energies between two Cs atoms and be- 
tween Nb and Cs are given in Table [Til Our total energy 
calculations suggest that pairing between Cs atoms is 
not energetically favorable. On the other hand Cs atoms 
can reduce their energy if they are adsorbed close to a 
Nb atom rather than on top of a Mo atom. This can 
be easily explained with the electrostatic attraction be- 
tween the two oppositely charged defects. This case of co- 
doping with Cs and Nb is particularly interesting, since 
it can be used as a strategy for fabricating p-n junctions, 
or more generally devices that require both n-type and 
p-type conductivity. Importantly the DOS of a super- 
cell where a Mo-substitutional Nb impurity is bonded to 
a Cs adatom appears essentially identical to that of an 
undoped M0S2 monolayer (see Fig. [s]), indicating that 
charge compensation is very effective. As a consequence 
our results suggest that an effective strategy for creating 
p-n hetero-j unctions must prevent the formation of Cs- 
Nb bonding by suppressing the migration of Cs on the 
surface. 

In summary we find that adsorbed alkali metals release 
their valance s electron to the M0S2 conduction band, 
and that the adsorption energies are large and negative, 
meaning that the doping will occur easily. Therefore, ad- 
sorbed alkali metals, with the exception of H, appear as 
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FIG. 9: Schematic diagrams of the electric double layer tran- 
sistor operation with an ionic liquid electrolyte. Here, "S" in- 
dicates the source, "D" indicates the drain, and "G" indicates 
the gate. The solid circles denote the ions in the liquid, where 
the different colors represent cations and anions, (a) When no 
gate voltage is applied (Vg=0), cations and anions are uni- 
formly distributed and both adsorbed at the interface of the 
semiconductor with equal probability, (b) When a finite gate 
bias is present (Vg = v), cations or anions adsorb predomi- 
nantly at the gate electrode, depending on the bias polarity, 
and the oppositely charged ions adsorb predominantly on the 
semiconductor. These adsorbed ions lead to an accumulation 
of a screening charge at the semiconductor surface, which im- 
plies a large surface carrier concentration. 



ideal candidates for doping n-type M0S2. To put this re- 
sult in context we remark that various experimental stud- 
ies have demonstrated the possibihty of dope n-type bulk 
M0S2 by intercalating alkali metala^. Here we demon- 
strate that the same is also possible at the single layer 
level. 
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FIG. 10: (Color online) The optimized geometries for (a, b) 
NHJ and (c, d) BF4 ions adsorbed on a M0S2 monolayer. In 
(a) and (c) the top view is shown, whereas the side view is in 
(c) and (d). Pink spheres indicate H, green spheres N, cyan 
spheres F and red spheres B. 




E - Ep (eV) 



Adsorption of molecular ions 

Electric double layer (EDL) transistors with ionic liq- 
uids (ILs) employed as gate dielectrics have recently 
emerged as promising devices, where the electrical prop- 
erties of a s olid can be controlled by electrostatic carrier 
dopin^25|49| A schematic diagram for such a device is 
shown in Fig. |9] By applying a gate bias to the ionic 
liquid, an electric double layer is formed between the in- 
terface of the liquid and the solid (see Fig. [9]) . Here the 
polarity of the molecular ions on the surface of the semi- 
conductor can be reversed by changing the polarity of the 
gate bias. The accumulation of high charge carrier densi- 
ties is possible ('^10"'^* cm~^) and this can be modulated 
by the applied gate voltage. In some case the accumu- 
lation can be so large to lead to phase transitions in the 
solid. Metal-superconductor and metal-insulator transi- 
tions have been demonstrated in KTiOj^ and V02^ , 
respectively. The operation of highly flexible M0S2 thin- 
film transistors in an ionic liquid gated dielectric has been 
demonstrated recently ^^. Here we discuss the effects of 
the adsorbed molecular ions on the M0S2 monolayer. 

A variety of different ionic liquids are experimentally 
available. Here we consider NH^ as cation and BF4 as 
anion, and we place them on the 5x5 M0S2 supercell 



FIG. 11: (Color online) DOS for (a) NH^ and (b) BFJ 
molecule adsorbed on M0S2. The results reported here are 
for a 5x5 supercell. The blue dashed line indicates the Fermi 
energy, the colored shaded areas indicate the DOS projected 
over the adsorbates. Note that the DOSs have the Fermi level 
at (Ep = 0) but they are aligned according to their valence 
and conduction bands. 



density of ^10^"' cm~^. We find that the most energet- 
ically favorable adsorption site is Tg for the cation and 
Tmo for the anion. In Fig. 11 a) the DOS is shown for 



(see Fig. 10 1. Such geometry corresponds to an impurity 



NH4 adsorption. It can be seen that the NH4 releases one 
electron into the M0S2 conduction band, resulting in a 
n-type M0S2 monolayer channel. In contrast, when BF4 
is adsorbed, one electron is transferred from M0S2 to the 
molecule [see Fig. ITTVb)]. As a consequence the Fermi 
energy shifts below the M0S2 VBM, so that the system 
becomes p-type. This also confirms the fact that the ad- 
sorbed molecules are indeed ionized to NH^ and BF^ 
when put on M0S2. In a liquid gated transistor arrange- 
ment, the anions/cations are brought close to the surface 
of the channel by applying a bias at the gate. Therefore, 
for this system, depending on the polarity of the gate 
bias the M0S2 is indeed predicted to have a n- or p-type 
conducting channel. Clearly this effect is maximized for 



a single layer of M0S2 when compared to thicker struc- 
tures. In fact, conductivity measurementa^Sl on M0S2 
show that EDL gating becomes less efficient when go- 
ing from a thin layer to bulk. Moreover, depending on 
the polarity of the bias, it is found that the conductivity 
indeed switches from n-type to p-type. 



Robustness of the results against the choice of 
XC functional:HSE06 



Finally, in order to verify that the impurity level align- 
ment presented in the previous sections is robust against 
the level of approximation taken for the DFT exchange- 
correlation energy, we here repeat our calculations by 
using the HSE06 functional. The M0S2 lattice constant 
obtained with HSE06 is 3.155 A, which is only slightly 
larger than the LSDA value of 3.137 A. For a pristine 
M0S2 monolayer at this relaxed lattice constant HSE06 
exhibits a direct band-gap of 2.2 eV (compared to a LDA 
band-gap of 1.86 eV). This is rather similar to the band- 
gap obtained by applying the atomic self-interaction cor- 
rection scheme^. 




FIG. 12: (Color online) DOSs for the doped monolayer M0S2, 
calculated with the HSE06 functional. We report results for 
(a) Nb substitutional at the Mo site, and for (b) Cs, (c) NH^ 
and (d) BFJ adsorbed on the M0S2 surface. The blue dashed 
line indicates the Fermi energy {Ep = 0), the colored shaded 
areas indicate the DOS projected on the adsorbates. In the 
plot we have aligned all the conduction and valence bands. 



In Fig. [12] we report our results for four representative 
doping strategies, namely for Nb substitutional doping, 
as well as for Cs, NH4 and BF4 adsorption. These corre- 
spond to the most promising dopants for either n-type or 



p-type doping as obtained at the LSDA level. The cal- 
culations are carried out with the same structures used 
in the previous sections, in order to make a direct com- 
parison between results obtained with the two function- 
als used. We note, however, that the M0S2 band-gap 
is very sensitive to changes in the lattice constanlp2l, so 
that the IISE06 band-gap at the LSDA lattice constant 
of 3.137 A becomes 2.31 eV. We hnd that for all the 
four cases the main difference between the LSDA and the 
HSE06 results is just a change in the band-gap, while the 
impurity levels with respect to the relative Fermi energy 
are placed at the same position. Therefore we can con- 
clude that the four most promising candidates for doping 
M0S2 monolayers remain those identified by our initial 
LSDA analysis. 



IV. CONCLUSION 

We have carried out a systematic study on the changes 
to the electronic structure of a M0S2 monolayer induced 
by various dopants. This study has the aim of identi- 
fying possible doping strategies for M0S2 monolayers to 
be used in devices manufacturing. We have considered 
first substitutional doping at both the Mo and the S site, 
as well as doping by adsorption. In general, S substitu- 
tion with non-metals and Mo substitution with transi- 
tion metals create deep donor levels inside the band-gap 
of the M0S2 monolayer for most of the dopants consid- 
ered. However, with substitutional doping we find that 
it is possible to obtain p-type M0S2 by replacing a Mo 
atom with Nb. In contrast n-type doping does not ap- 
pear possible since we have found deep donor levels for 
all the substitutions investigated (either at the Mo or the 
S site). More promising n-type doping can be achieved 
by adsorbing alkali metals on the surface of M0S2. Fi- 
nally, as a last class of dopants, we have considered the 
adsorption of ionic molecules, which occurs during the 
electric double layer formation in M0S2 embedded in an 
ionic liquid. These show high potential for inducing large 
carrier concentrations within the M0S2 monolayer, in the 
form of both electrons or holes. 
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